Quantum state engineering, purification, and number resolved photon detection with 

high finesse optical cavities 
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We propose and analyze a multi-functional setup consisting of high finesse optical cavities, beam 
splitters, and phase shifters. The basic scheme projects arbitrary photonic two-mode input states 

onto the subspace spanned by the product of Fock states |n)|n) with n = 0, 1, 2, This protocol 

does not only provide the possibility to conditionally generate highly entangled photon number states 
as resource for quantum information protocols but also allows one to test and hence purify this type 
of quantum states in a communication scenario, which is of great practical importance. The scheme 
is especially attractive as a generalization to many modes allows for distribution and purification 
of entanglement in networks. In an alternative working mode, the setup allows of quantum non 
demolition number resolved photodetection in the optical domain. 

PACS numbers: 42.50.Dv, 42.50.Pq, 03.67.-a 



I. INTRODUCTION 

Light plays an essential role in quantum communica- 
tion and is indispensable in most practical applications, 
for example quantum cryptography. Photons are attrac- 
tive carriers of quantum information because the interac- 
tions of light with the surroundings are normally weak, 
but for the same reason it is generally difficult to pre- 
pare, manipulate, and measure quantum states of light 
in a nondestructive way. Repeated interactions provide 
a method to increase the effective coupling strength be- 
tween light and matter, and the backreflection of light in 
a cavity thus constitutes an interesting tool, in particu- 
lar, because experiments are currently moving into the 
strong coupling regime [lj-Q , where coherent dynamics 
takes place on a faster time scale than dissipative dy- 
namics. 

In this paper we propose a versatile setup consisting of 
an array of cavities and passive optical elements (beam 
splitters and phase shifters), which allows for quantum 
state engineering, quantum state purification, and non- 
destructive number resolving photon detection. The 
setup builds on two basic ingredients: The Hong-Ou- 
Mandel interference effect [f| generalized to input pulses 
containing an arbitrary number of photons and the pos- 
sibility of projection onto the subspace of even or the 
subspace of odd photon-number states by use of cavity 
quantum electrodynamics in the strong coupling regime. 

Regarding quantum state engineering, the basic setup 
provides a possibility to conditionally generate photon- 
number correlated states. More specifically, the setup 
allows us to project an arbitrary photonic two-mode in- 
put state onto the subspace spanned by the state vectors 

\n)\n) with n = 0,1,2, We denote this subspace by 

S. The scheme is probabilistic as it is conditioned on a 
specific measurement outcome. The success probability 
equals the norm of the projection of the input state onto 
S and is thus unity if the input state already lies in S. 



In other words, the setup may be viewed as a filter |6j, 
which removes all undesired components of the quantum 
state but leaves the desired components unchanged. We 
may, for example, use two independent coherent states 
as input and obtain a photon-number correlated state as 
output. 

Photon-number correlated states, for example 
Einstein-Podolsky-Rosen (EPR) entangled states [3], 
are an important resource for quantum teleportation 
entanglement swapping [I3 M l 5| . quantum key 
distribution [r| [TtJ , and Bell tests [l8l Il9l] . In practice, 
however, the applicability of these states is hampered by 
noise effects such as photon losses. Real- world applica- 
tions require therefore entanglement purification. The 
proposed setup is very attractive for detection of losses 
and can in particular be used to purify photon-number 
entangled states on site. If a photon-number correlated 
state, for example an EPR state, is used as input, the 
desired state passes the setup with a certificate, while 
states which suffered from photon losses are detected 
and can be rejected. 

Photon losses are an especially serious problem in 
quantum communication over long distances. It is not 
only a very common source of decoherence which is hard 
to avoid, but also typically hard to overcome. The on- 
site purification protocol mentioned above can easily be 
adopted to a communication scenario such that it allows 
for the purification of a photon-number correlated state 
after transmission to two distant parties. 

Purification of two mode entangled states has been 
shown experimentally for qubits [2(1 [2l| and in the con- 
tinuous variable (CV) regime [H, HH • (CV- entanglement 
purification is especially challenging |2#j2a |. Neverthe- 
less, several proposals have been made to accomplish this 
task [27l - l32| , and very recently Takahashi et al. succeeded 
in an experimental demonstration [33[.) A special advan- 
tage of our scheme lies in the fact that it does not only 
allow for detection of arbitrary photon losses, but is also 
applicable to many modes such that entanglement can 



be distributed and purified in a network. 

With a small modification, the basic setup can be used 
for number resolved photon detection. The ability to 
detect photons in a number resolved fashion is highly 
desirable in the fields of quantum computing and quan- 
tum communication. For example, linear optics quantum 
computation relies crucially on photon number resolving 
detectors [3J-|36j. Moreover, the possibility to distin- 
guish different photon-number states allows for condi- 
tional state preparation of nonclassical quantum states 
[37M39I ] . and plays a role in Bell experiments lioll and 
the security in quantum cryptographic schemes [4ll . |42| . 
Other applications include interferometr v 1431 a nd the 
characterization of quantum light sources [44j, |45|. 

Existing technologies for photon counting [46l459j such 
as avalanche photodiodes, cryogenic devices, and quan- 
tum dots typically have scalability problems and cannot 
reliably distinguish high photon numbers, destroy the 
quantum state of light in the detection process, or do 
not work for optical photons. Here, we present a non- 
destructive number resolving photo detection scheme in 
the optical regime. This quantum-non-demolition mea- 
surement of the photon number allows for subsequent use 
of the measured quantum state of light. An advantage of 
the counting device put forward in this work compared 
to other theoreti cal prop osals for QND measurements of 
photon numbers [60-65] is the ability to detect arbitrar- 
ily high photon numbers with arbitrary resolution. The 
scheme is based on testing successively all possible prime 
factors and powers of primes and the resources needed 
therefore scale moderately with (width and mean of the) 
photon number distribution. In particular, a very precise 
photon number measurement can be made even for very 
high photon numbers by testing only few factors if the 
approximate photon number is known. 

The paper is structured as follows. We start with a 
brief overview of the main results in Sec. |n] In Sec. IIII| 
we explain how the conditional projection onto S can 
be achieved and discuss some properties of the proposed 
setup in the ideal limit, where the atom-cavity coupling is 
infinitely strong and the input pulses are infinitely long. 
In Sec. IIV1 we show that a modified version of the setup 
can act as a non-destructive photon number resolving 
detector, and in Sec. |V] we investigate the possibility to 
use the setup to detect, and thereby filter out, losses. In 
Sec. IVI[ we consider the significance of finite input pulse 
length and finite coupling strength, and we obtain a sim- 
ple analytical expression for the optimal choice of input 
mode function for coherent state input fields. Section IVIII 
concludes the paper. 



II. OVERVIEW AND MAIN RESULTS 

The most important ingredient of the proposed setup 
is the possibility to use the internal state of a single atom 
to control whether the phase of a light field is changed by 
7r or not [66j . The basic mechanism, which is explained in 
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FIG. 1: (Color online) A single atom with level structure 
as shown in (a) is placed in a cavity in the strong coupling 
regime. The light field, which is on resonance with the cavity, 
couples the ground state level \f) resonantly to the excited 
state |e), and the state |e) decays to the state |t) through 
spontaneous emission at a rate F. (b) If the atom is in the 
state |4_) , the incident field is not affected by the presence of 
the atom, and for a sufficiently slowly varying input pulse, 
the interaction with the resonant cavity changes the phase of 
the light field by it. (c) If the atom is initially in the state 
|t), the possibility of spontaneous emission prevents the light 
field from building up inside the cavity (provided the photon 
flux of the input beam is not too high), and the incoming field 
is reflected from the input mirror without acquiring a phase 
shift. This transformation is insensitive to the precise values 
of g, F, and the cavity decay rate as long as the system is in 
the strong coupling and weak driving regime. 



Fig- H] has several possible applications, including prepa- 
ration of superpositions of coherent states [67j . continu- 
ous two-qubit parity measurements in a cavity quantum 
electrodynamics network (68j . and low energy switches 
[69| . Concerning the experimental realization, basic in- 
gredients of the scheme such as trapping of a single atom 
in a strongly coupled cavity and preparing of the initial 
atomic state have been demonstrated experimentally in 
[73], where the decrease in cavity field intensity for an 
atom in the state |f) compared to the case of an atom 
in the state |t) is used to subsequently measure the state 
of the atom. State preparation and readout for a single 
atom in a cavity have also been demonstrated in [71j . 
Another promising candidate for an experimental real- 
ization is circuit quantum electrodynamics. See, for in- 
stance, (72J for a review. 

The generation of quantum superposition states can be 
achieved as follows. The atom is initially prepared in the 
state (|t) + \l))/V2, and the input field is chosen to be 
a coherent state \a). After the interaction, the combined 
state of the atom and the light field is proportional to 
|a)|t) +\-a)\i) oc (\a) + | - a))(|t) + It)) + (|a) - | - 
a)) (|t) — |t)), and a measurement of the atomic state in 
the basis |±) = (|t) ± \i))/V2 projects the state of the 
light field onto the even \a) + \ — a) or the odd \a) — | — a) 
superposition state. More generally, the input state 
yt c n |n), where |n) is an n-photon Fock state, is trans- 
formed into the output state 5Z„(l/2 ± (— l)"/2)c„|n), 
i.e., the input state is conditionally projected onto either 
the subspace spanned by all even photon-number states 
or the subspace spanned by all odd photon-number states 
without destroying the state. 
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With this tool at hand, we can project an arbi- 
trary two-mode input state onto the subspace S — 

span(|n)|n)), n = 0, 1,2, If two modes interfere at a 

50:50 beam splitter, a state of form |n)|n) is transformed 
into a superposition of products of even photon-number 
states. If we apply a 50:50 beam splitter operation to the 
input state, project both of the resulting modes condi- 
tionally on the subspace of even photon-number states, 
and apply a second 50:50 beam splitter operation, the 
input state is thus unchanged if it already lies in S, but 
most other states will not pass the measurement test. 
To remove the final unwanted components, we apply op- 
posite phase shifts to the two modes (which again leaves 
\n)\n) unchanged) and repeat the procedure (as shown in 
Fig. [2]). For an appropriate choice of phase shifts, the de- 
sired state is obtained after infinitely many repetitions. 
In practice, however, a quite small number is typically 
sufficient. If, for instance, the input state is a product of 
two coherent states |a)|a) with \a\ 2 = 4, the fidelity of 
the projection is 0.573 for one unit, 0.962 for two units, 
and 0.999998 for three units. The scheme is easily gen- 
eralized to an M mode input state. In this case, we first 
project modes 1 and 2 on S, modes 3 and 4 on S, etc, 
and then project modes 2 and 3 on S, modes 4 and 5 on 
S, etc. 

The setup can also be used as a device for photon num- 
ber resolving measurements if the phases applied between 
the light-cavity interactions are chosen according to the 
new task. Each photon- number state \n) sent through the 
array leads to a characteristic pattern of atomic states. 
As explained in section llV Al one can determine the pho- 
ton number of an unknown state by testing the prime 
factors and powers of primes in the range of interest in 
subsequent parts of the array. The scheme scales thereby 
moderately in the resources. Three cavity pairs suffice 
for example for detecting any state which is not a mul- 
tiple of three with a probability of 93.75%. However, in 
this basic version of the counting scheme, the tested pho- 
ton state may leave each port of the last beam splitter 
with equal probability. Deterministic emission of the un- 
changed quantum state of light into a single spatial mode 
is rendered possible if we allow atoms in different cavi- 
ties to be entangled before the interaction with the field 
(see section HVB[) . More generally, the proposed scheme 
allows to determine the difference in photon numbers of 
two input beams without changing the photonic state. 

The correlations in photon number between the two 
modes of states in S facilitate an interesting possibility 
to detect photon losses. To this end the state is projected 
onto S sl second time. If photon loss has occurred, the 
state is most likely orthogonal to 5, in which case we 
obtain a measurement outcome, which is not the one we 
require in order to accept the projection as successful. On 
the other hand, if photon loss has not occurred, we are 
sure to get the desired measurement outcome. We note 
that loss of a single photon can always be detected by this 
method, and the state can thus be conditionally recov- 
ered with almost perfect fidelity if the loss is sufficiently 



small. We can improve the robustness even further, if 
we use an M-mode state. It is then possible to detect 
all losses of up to M — 1 photons, and even though it 
is M times more likely to lose one photon, the probabil- 
ity to lose one photon from each mode is approximately 
(Mp) M , where p is the probability to lose one photon 
from one mode and we assume Mp < 1. In a situation 
where many photon losses are to be expected, this pro- 
cedure allows one to obtain photon-number correlated 
states with high fidelity, although with small probability. 

We can also distribute the modes of a photon-number 
correlated state to distant parties, while still checking for 
loss, provided we send at least two modes to each party. 
As the proposed scheme can be used as a filter prior to the 
actual protocol it has an important advantage compared 
to postselective schemes. If the tested entangled state 
is for example intended to be used for teleportation, the 
state to be teleported is not destroyed in the course of 
testing the photon-number correlated resource state. 

The dynamics in Fig. [T] requires strong coupling, a suf- 
ficiently slowly varying mode function of the input field, 
and a sufficiently low flux of photons. To quantify these 
requirements, we provide a full multi-mode description of 
the interaction of the light with the cavity for the case of 
a coherent state input field in the last part of the paper. 
We find that the single atom cooperativity parameter 
should be much larger than unity, the mode function of 
the input field should be long compared to the inverse of 
the decay rate of the cavity, and the flux of photons in 
the input beam should not significantly exceed the rate 
of spontaneous emission events from an atom having an 
average probability of one half to be in the excited state. 
We also derive the optimal shape of the mode function 
of the input field (Eq. (|27| ). when the mode function is 
only allowed to be nonzero in a finite time interval. 



III. NONDESTRUCTIVE PROJECTION ONTO 
PHOTON-NUMBER CORRELATED STATES 

The proposed setup for projection of an arbitrary two- 
mode input state onto S is sketched in Fig. [5] We denote 
the field annihilation operators of the two input modes 
by a and b, respectively. The total transformation corre- 
sponding to one of the units consisting of a beam splitter, 
a set of cavities, and a second beam splitter, conditioned 
on both atoms being measured in the state |+) after the 
interaction, is given by the operator [/' PU, where 



U = exp 



and 



P=Y2J2 \ 2n )( 2n \ ® |2m)(2r 



(1) 



(2) 



n— m— 



As explained above, the Hong-Ou-Mandel effect ensures 
that U'PU\n)\n) = \n)\n), while most other possible 
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FIG. 2: (Color online) The first three units of the proposed setup to conditionally project an arbitrary two-mode input state 
onto the subspace spanned by the state vectors |n)|n), n = 0, 1, 2, . . .. All atoms are prepared in the state (|t) + i))/\/2 before 
the interaction with the field, and the desired projection occurs in the limit of infinitely many units conditioned on all atoms 
being in the state (|t) + IJ,})/"^ after the interaction. As explained in the text, a small number of units will typically suffice 
in practice. For later reference we label the beam splitters as BSi and the cavities as Ci. 



components of the input state are removed through the 
conditioning, for instance all components |^)|to) with 
n + to odd. There are, however, a few exceptions, 
since all states of the form U'\2n)\2m) , n = 0,1, 2,..., 
to = 0, 1, 2, . . ., are accepted. The phase shifts between 
the U'PU units are represented by the operator 



U,), = exp i(f> (at a 



ftb 



(3) 



which leaves states of the form unchanged, while 

states of the form |n)|m) with n ^ to acquire a phase 
shift. 

For a setup containing N+l units, the complete condi- 
tional transformation is thus represented by the operator 

6 N = UiPUUfrrfPU ■■■UfctfPUU^UiPU (4) 



A' 



= L/ t PC/]Jcos[^(a t a-fe t 6)], 



(5) 



not measured before the interaction with the light field 
is completed. The setup using an array of cavities as in 
Fig. [5] may thus be very compact even though the pulses 
are required to be long. (Note that it would also be pos- 
sible to use a single pair of cavities and atoms repeatedly 
in a fold-on type of experiment; however in that case the 
compactness would be lost due to the need for long de- 
lay lines necessary to measure and re-prepare the atoms 
before they are reused.) 

A natural question is how one should optimally choose 
the angles <pi to approximately achieve the projection 
with a small number of units. To this end we define 
the fidelity of the projection 



Fn — 



l(<M</WI 5 



(■>pN\'>fiN)(lp 00 \4> 



(V'oolV'oo) 



'N 



(6) 



where W PU in the last line commutes with the product 
of cosines. For N — > oo, the product of cosines van- 
ishes for all components of the input state with different 
numbers of photons in the two modes if, for instance, all 
the (pi's are chosen as an irrational number times 7T. We 
note that even though we here apply the two-mode oper- 
ators one after the other to the input state corresponding 
to successive interactions of the light with the different 
components of the setup, the result is exactly the same 
if the input pulses are longer than the distance between 
the components such that different parts of the pulses in- 
teract with different components at the same time. The 
only important point is that the state of the atoms is 



as the overlap between the unnormalized output state 
IV'w) = Ojv|^in) after N + l units and the projection 
\ipoo) of the input state \ipia) onto the subspace S. The 
last equality follows from the fact that IV'jv) = IV'oo) + 
\i/>±), where \ip±) lies in the orthogonal complement of S. 
Maximizing F N for a given |^ ln ) = E„ H m c ™l n )l m ) 
thus corresponds to minimizing 



oo oo N 

(iPn\iPn) = X! X! c «m JJcos 2 [0i(n - to)] 



n=0 m—0 2—1 



x (ilto.\utpu\n)\m), (7) 
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i.e., we would like to find the optimal solution of 



d(ip N \ipN) 



^ c„ m sin[20j(n - m)](n - m) 



n— m— 



JV 



JJcos^iin - m)]{>ipin\rf PU\n)\m) = 0. (8) 



A set of solutions valid for any input state can be 
obtained by requiring sin[20j(n — m)} Ylntj cos 2 [<pi(n — 
to)] = for all even values of n — m (note that 
W PU\n)\m) = for n + m odd). Within this set 
the optimal solution is 4>j = 2~ J x tt/2. It is inter- 
esting to note that choosing one of the angles to be 
2~ J x 7r/2, j e {1,2,..., N}, all terms with n - m = 
2 J x (±1, ±3, ±5, . . .) are removed from the input state 
according to ([5]). When all angles are chosen according 
to 4>j = 2~- 7 x 7r/2, it follows that \4>n) only contains 
terms with n — to = q2^ N+1 \ q = 0, ±1, ±2, . . ., which 
may be a useful property in practical applications of the 
scheme. Even though this is not necessarily optimal with 
respect to maximizing F for a particular choice of input 
state, we thus use the angles <f>j — 2~ J x tt/2 in the fol- 
lowing, except for one important point: If the input state 
satisfies the symmetry relations c nm = c mn , it turns out 
that the operator U'PU by itself removes all terms with 
n — m = ±2, ±6, ±10, . . ., i.e., we can choose the angles 
as (f>j = 2~ 3 x 7r/4, and \i/jn) only contains terms with 
n-m = q2 ( - N+2 \ q = 0,±1,±2,.... For N = 2, for 
instance, only terms with n — to = 0, ±16, ±32, . . . con- 
tribute. 

In Fig. [31 we have chosen the input state to be a prod- 
uct of two coherent states with amplitude a and plotted 
the fidelity JB} as a function of \a\ 2 for different num- 
bers of units of the setup. Even for \a\ 2 as large as 10, 
the fidelity is still as high as 0.9961 for N = 2, and the 
required number of units is thus quite small in practice. 
The figure also shows the success probability 



Pn = (tPn\4>n) 



(9) 



for N — > oo. For \a\ 2 = 10, for instance, one should 
repeat the experiment about 11 times on average before 
the desired measurement outcome is observed. 



IV. PHOTON NUMBER RESOLVING 
MEASUREMENT 

In this section, we show how a photon number mea- 
surement can be implemented using a modified version 
of the setup introduced in the previous section. The key 
idea is explained in subsection IIV Al where we describe 
the basic photo-counting scheme. In subsection IIVB1 
this protocol is extended to allow for a QND measure- 
ment of photon numbers. 




FIG. 3: (Color online) Fidelity (Eq. (JBJ) as a function of the 
expectation value of the number of photons in one of the input 
modes for \tpi n ) = |q)|o!} and setups with one, two, three, 
and infinitely many units. The angles are chosen as cj>j = 
2~i x 7r/4. The dotted line labeled is the probability in the 
limit of infinitely many units to actually obtain the required 
measurement outcome, i.e., all atoms in (|t) + \\,))/\/~2 after 
the interaction with the field. 



A. Number resolving detection scheme 

In the following, we analyze the setup shown in Fig. [2] 
when the input is a product of an n-photon Fock state in 
the lower input beam and a vacuum state in the upper 
input beam. Since the setup contains a series of beam 
splitters, it will be useful to define A = (a> — w) / \/2 and 
B = (fit + St)/V2, such that a T |0) -> i|0) and S f |0) -> 
B\0) at beam splitters BSi, BS 3 , BS 5 , . . ., and A\0) -> 
ot|0) and B\0) -> S+|0) at beam splitters BS 2 , BS 4 , BS 6 , 

As before, all atoms are initially prepared in the state 
|+) and will after the interaction with the field be mea- 
sured in the |±) basis. When we start with an n photon 
state, there are only two possible outcomes of the mea- 
surement of the atoms in the cavities labeled Ci and C2 
in Fig. [5] dependent on whether n is even or odd. In the 

even case, the two atoms can only be in | + +) or | ) 

and in the odd case | — h) or | H — ). To handle the odd 
and even case at the same time we denote | + +) , | — h) 

as \B + ) and | ), | H — ) as \B-). A measurement of 

\B + ) indicates an even number of photons in the &-beam, 
resp., \B—) an odd number of photons in the &-beam. 

We start with the state |n)|0) = ^(a T )™|0)|0) as in- 
put. After the beam splitter BSi, the state has changed 
and interacts with the two atoms 



into 



in the cavities Ci and C2. By measuring the atoms 
in the |±) basis, the state is projected into the sub- 
space of an even or odd number of photons in the b 
path. The photon state after the measurement can be 
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written as |6j 



-7= (A" ± B n )|0)|0), where is the state with 

an even (odd) number of b photons and corresponds to 
the measurement result \B+) Note that this first 

measurement result is completely random. 



After BS2 the state simplifies to 



[(at)« ± 



(&t)™]|0)|0). Now due to the phase shifters the two modes 
pick up a relative phase of 2<pin so that the state is given 



by 



,i<f>in 



(at)" ± 



-i<fi\n 



/2 n\ 

BS3 we have the state 
which is equal to 1 

1 ^ c i<t>-in 



(fet)»]|0)|0). 



Finally, after 

_l ^ e i<Pm^n ± e -#i"Bn)|0)|0), 

n ,(e* n ±e-^ in )(i"+B")|0)|0) + 
*")(> - B n )|0)|0). This can also 
(e^ 1 " ± e-^ in )/2|6 + ) + (e^ 1 " 



2%/2 n! 

be rewritten as (e" f ' 1 " ± e~' vl ")/2|6 + ) + (e"'' 1 " =p 
e~** >1 ")/2|&_). So the result of measuring the state 
of the atoms in cavities C3 and C4 will be \B±) with 
probability p + = cos(0in) 2 and \B^) with probability 
p_ = sin(0in) 2 . Since the state is again projected into 
one of the two states \b±) we can repeat exactly the same 
calculations for all following steps. 

Whereas the first measurement result was completely 
random, all following measurement results depend on n 
and the previous measurement outcome, i.e., with proba- 



bility pi = cos( 



the (i + l)-th measurement result is 



the same as the i-th result, and with probability sin(<fon) 
the measurement result changes and the state changes 
from \b±) to |6=p). If the number of units is infinite (or 
sufficiently large) and we have chosen all phases equal 
as <j), then we can guess from the relative frequency with 
which the measurement result has switched between \B+) 
and \B—) the number of photons with arbitrary precision 
for all photon numbers n < ?. n~ arccos( v / /)/0, where 

/ = N same /(N samc + TV d 

iffcrcnt 

) and iVsamc (TV di 

ffcrcnt 

the number of cases, where the measurement outcome is 
the same (not the same) as the previous measurement 
outcome. 

Measuring this relative frequency with a fixed small 
phase is not the optimal way to get the photon number. 
We propose instead the following. Let us use a setup 
with a total of TV + 1 units and choose the phases to be 
(j>i = 2 l ~ 1 ir/no, i = 1, 2, . . . , N, for an arbitrarily chosen 
value of no 6 N and let us calculate the probability p(n) 
that the measurement results are all the same, 



A' 



N-l 



p( n ) = Yi Pi = n c ° s ( 

»=1 i=0 ^ n ° 



(10) 



This probability is equal to one for all photon numbers 
that are a multiple of hq and goes to zero otherwise in 
the limit of infinitely many units of the setup. This way 
we can measure whether the photon number is a multiple 
of no- 

For example, for no = 3 and N + 1 = 3 we detect any 
state which is not a multiple of three with a probability 
of at least q = 93.75%, resp., q = 99.61% for N + 1 = 5, 




FIG. 4: (Color online) Probabilities p(n) (Eq. flT0j) for n = 
20 and N + 1 = 5 resulting in q = 94.49% (Eq. JTT)). For 
TV + 1 = 7, q = 99.66%. 




FIG. 5: (Color online) Probabilities p(n) for no — 100 and 
N + 1 = 8 resulting in q = 96.33%. For N + 1 = 10, q = 
99.95%. 



where 



1 — max p(n). 

n^0,no,2no,... 



(11) 



For no — 4, already N + 1 = 3 is sufficient to achieve 
q = 100%. For 7j = 5 and N + 1 = 3 we have q = 
93.75%, which increases to q = 99.61% for TV + 1 = 5. 
In Fig. U and [S] we have shown p(n) for no = 20 and 
no = 100. The number TV needed to get a good result 
typically scales logarithmical with no, e.g. for no = 1000 
already TV + 1 = 11 is enough to reach q = 99.95%. 

Given an unknown state we can test all possible prime 
factors and powers of primes to identify the exact pho- 
ton number. If, e.g., we have a state consisting of 
to 10 photons the following factors have to be tested 
n = 2, 3, 4, 5, 7, 8, 9 (where 2 does not need to be checked 
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separately) . 24 measurement results are sufficient to test 
all factors with a probability over 99%. 

If we require reliable photon number counting for n 
ranging from to n max , for large n max , all primes and 
power of primes that are smaller than n max need to be 
tested. This number can be bounded from above by 
"■max- All u max tests are required to work with high 
probability. To this end each single test needs to suc- 
ceed with a probability better than q > 1 — l/n max . 
It can be checked numerically that this is the case if 
N = 21og(n max ), leading to a photon counting device 
with reliable photon detection up to n max using an array 
consisting of less than 2n max log(n max ) basic units. 

Note that this setup does not destroy the photonic in- 
put state but changes |n)|0) randomly into -^[|*n)|0) ± 
1 0) | n.) ] , i.e., the photons leave the setup in a superposition 
of all photons taking either the a-beam or the &-beam. 
This can not be changed back into |n)|0) by means of 
passive optical elements. The output state - a so-called 
./V00-/V state - is, however, a very valuable resource for 
applications in quantum information protocols and quan- 
tum metrology [73|- 



B. Non-destructive number resolving detection 
scheme 



For a non-demolition version of the photon number 
measurement we use the basic building block depicted in 
Fig. |5J The atoms in the two upper cavities are initially 
prepared in an entangled state \cj) + ) = (|tt) + |44))/v2) 
and the atoms in the lower cavities are also prepared 
in the state | </>+). This can, for instance, be achieved via 
the parity measurement scheme suggested in [68| . During 
the interaction the upper and the lower atoms will stay 
in the subspace spanned by \<f>±) = (|tt}=M-U-})/'V / 2. The 
state changes between and \4>-) each time one of the 
two entangled cavities interacts with an odd number of 
photons. As in the previous subsection we define \B±) to 
handle the even and the odd case at the same time. For 
n even, \B+) = \cf)+)\(j)+) and \B~) = \<j)-)\(f)-), while 
in the odd case we define \B + ) — |</>_)|</> + ) and \B-) = 

\<P+)\4>-)- 

Using the same notation as before, the state is trans- 
formed as follows, when we go through the setup from 
left to right. The initial state is 



i^ 1 ) = 4=(« t rio)io)i0 + )i0+), 

V n\ 



and after the first beam splitter we have the state 



\i> 2 ) = -Li"|o>|o)|0 + >I4>+) = i(IM+l&-))l<MI^ 



The interaction with the first two cavities leads to the 



IV>1> IV'2> 

Ci 



\ \<t>±) 



\4>r) \4>») 




FIG. 6: (Color online) Basic building block of the non- 
demolition photon number resolving detection scheme. The 
dashed circles and the wavy lines indicate the entanglement 
between the atoms in cavities Ci and C3 and between the 
atoms in cavities C2 and C4. The gray box is either a cavity 
or a mirror. (In the latter case, we may as well send the light 
directly from the cavities C3 and C4 to the first two cavities 
of the next block.) 



state 



\fr) = ±\b + )\B + ) + ±=\b-)\B-) 



1 



2-Jn\ 



(A n + B n )\0)\0)\B^ 



1 



2s/n\ 



(A n - B n )\0)\0)\B_ 



The second beam splitter transforms the state into 

|V, 4 ) = -l={[(at)" + (St)"]|0)|0>|S + ) 
2Vn! 



+ [(at)™-(6t)»]| )|0)|S_)}. 

The two modes pick up a relative phase shift of 2<^n at 
the phase shifters so that 

|V 5 > = -^={[e^(atr + e-*«(St)"]|0)|0)|B + ) 
2\m! 

After the third beam splitter we get 

|V> 6 ) = ^^[(e^"i n + e-^"B")|0)|0)|B + ) 
2\/n\ 

+ ^n^n _ e -^"B Tl )|0)|0)|B_)] 
= i cos(<t>in)\b+)\B + ) + -J= sin(^n)|6_)|B + ) 

i 1 
+ — sin(^7i)|b+)|B_) + — cos(0 4 n)|&_)|S_). 

Note now that independent of whether n is even or odd, 
the interaction with the last two cavities turns the states 
\b±)\B±) into \b±) \<j>+) \<j>+) and the states \b±)\Bzp) into 
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\b±)\(j> -)[(/>-). The state is thus changed to 



|V>7) = ^cos(<M(IM + IO)l<MI<M 

= -±=A n \0)\0) [cos(4n)|0 + >|0+> +ism(^n)|0_>|0_) 
v n! 



and after the last beam splitter we have 



IVs 



7=(a T )"|0)|0) 

/n 

<g> [cos(^n) | (/>+)|(/> + ) +isux((f>in)\<p-)\<f)-)] . 



Note that the photonic modes are now decoupled from 
the state of the atoms. The photons will continue af- 
ter the final beam splitter unchanged in |n)|0) while the 
atoms still contain some information about the photon 
number. 

Note that \4>+) = (|++} + |--))/V2and |0_) = (| + 
— ) + 1 — h))/ V%- By measuring all atoms in the |±) basis 
we can easily distinguish between \4>±) by the parity of 
the measurements. The probability to obtain |(/> + )|0 + ) is 
cos(0n) 2 , and the probability to get \(j>J)\(j)J) is sin(<^n) 2 . 
After the measurement, all photons are found in the a 
beam for both outcomes. The probability for measuring 
| (/>-!_) | </>_i_) is the same as the changing probability in the 
previous setup such that we can do the same with a chain 
of the demolition free block. If we prefer also to get the 
parity information in each step, we can add an additional 
cavity at the end of each block as shown in Fig. [51 

More generally, the demolition free element leaves pho- 
tonic input states \ip) = ^t g i (" t )"(^) <? |0)|0)i where n 
photons enter through the lower and q photons enter 
through the upper port, unchanged. A calculation anal- 
ogous to the previous one shows that one obtains the 
atomic states |</>+)|</>+) and \4>J)\4>J) with probabilities 
cos(4>i(n—q)) 2 and sin(0i(n— g)) 2 respectively. This way, 
we can test for photon number differences n — q in two 
input states in the same fashion as for photon numbers 
in a single input beam described above. Similarly, one 
can project two coherent input states onto gener- 

alized photon-number correlated states ^2 n c n \n)\n — d) 
with fixed photon number difference d = 0, 1, 2.... 

In a realistic scenario we may be faced with photon 
losses. Both setups have a built-in possibility to detect 
loss of one photon. In the hrst case we get the parity of 
the total number of photons in every single measurement 
of a pair of atoms. If this parity changes at some place 
in the chain then we know that we lost at least one pho- 
ton. In the demolition free setup the valid measurement 
results are restricted to |^+)|^+) and | ) | <fr- ) . If we 
measure |0_)|</> + ) or \(f> + )\<j)-) we know that we lost a 
photon in between the two pairs of cavities. In addition, 
the optional cavity at the end of each block provides an 
extra check for photon loss. 



■g 0.6 - 




FIG. 7: (Color online) Projecting the input state \ipi n ) = 
| a) | a) onto the subspace S followed by a fractional loss R in 
both modes and a second projection onto S, the figure shows 
the fidelity between the states obtained after the first and 
the second projection onto S, the probability that the second 
projection is successful given that the first is successful, and 
the purity of the state after the second projection for two 
different values of la! 2 . 



V. FILTERING OUT LOSSES 

We next investigate the possibility to use a compari- 
son of the number of photons in the two modes to de- 
tect a loss. As in Sec. IIII1 we start with the input state 
iV'in) = |o:)|a) and use the proposed setup to project 
it onto the subspace S. We then use two beam split- 
ters with reflectivity R to model a fractional loss in both 
modes. After tracing out the reflected field, we finally 
use the proposed setup once more to project the state 
onto S. In Fig. [7J we plot the fidelity between the state 
obtained after the first projection onto S and the state 
obtained after the second projection onto S, the proba- 
bility that the second projection is successful given that 
the first is successful, and the purity of the state after 
the second projection. The second projection is seen to 
recover the state obtained after the first projection with 
a fidelity close to unity even for losses of a few percent. 
This is the case because a loss of only one photon will 
always lead to a failure of the second conditional projec- 
tion. The main contribution to the fidelity decrease for 
small R is thus a simultaneous loss of one photon from 
both modes. It is also interesting to note that the final 
state is actually pure for all values of R, which is a con- 
sequence of the particular choice of input state. Finally, 
we note that a single unit is sufficient to detect loss of a 
single photon, and for small R we thus only need to use 
one unit for the second projection in practice. 

Let us also consider a four mode input state IV'in) = 
|a)|a)|a)|Q;). As before we use the setup to project this 
state onto the subspace spanned by the state vectors 
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the following. For simplicity, we assume the input field 
to be a continuous coherent state 174] 



|{a ta (t)}> = exp J a in (i)a t (t)dt - J a* n (t)a{t)dt |0) 

(12) 

with mode function f- m (t) — a- m (t)/a, where \a\ 2 = 
J \ai n (t)\ 2 dt is the expectation value of the total number 
of photons in the input beam. The light-atom interaction 
is governed by the Hamiltonian 



H = ftg(cV + CT T c), tr:=|t)<e| 



and the decay term 



cr up 



(13) 



(14) 



FIG. 8: (Color online) Projecting the input state \ipi-n) = 
| a) | a) | a) | a) onto the subspace spanned by the vectors 
|n)|n}|n}|n), n = 0,1,2,..., followed by a fractional loss R 
in all four modes and a projection of modes 1 and 2 onto S 
and of modes 3 and 4 onto 5, the figure shows the fidelity 
between the states obtained after the first and the second 
set of projections, the probability that the projections onto S 
are successful given that the first projection is successful, and 
the purity of the state after the projections onto S for two 
different values of la! 2 . 



|n)|n)|n)|n), n = 0,1,2,.... We then imagine a frac- 
tional loss of R to take place in all modes. If two of 
the modes are on their way to Alice and the two other 
modes are on their way to Bob, we can only try to re- 
cover the original projection by projecting the former two 
modes onto 5 and the latter two modes onto S. The re- 
sults are shown in Fig. and again the curves showing 
the fidelity and the purity are seen to be very flat and 
close to unity for small losses. This scheme allows one to 
distribute entanglement with high fidelity, but reduced 
success probability. 



VI. DEVIATIONS FROM IDEAL BEHAVIOR 

So far we have considered the ideal limit of infinitely 
long pulses and an infinitely strong coupling. In this sec- 
tion, we use a more detailed model of the interaction of 
the light field with an atom in a cavity to investigate 
how long the input pulses need to be and how large the 
single atom cooperativity parameter should be to approx- 
imately achieve this limit. 



A. Optimal input mode function 

The backreflection of light in a cavity leads to a state 
dependent distortion of the shape of the mode function 
of the input field, and to study this effect in more detail 
we concentrate on a single cavity as shown in Fig. [1] in 



where g is the light-atom coupling strength, c is the 
annihilation operator of the cavity field, T is the de- 
cay rate of the excited state of the atom due to spon- 
taneous emission, and p is the density operator repre- 
senting the state of the atom and the cavity field. For 
g 2 Tr(c^c(t\p\t)) < (r/2) 2 , where Tr denotes the trace, 
the population in the excited state of the atom is very 
small, and we may eliminate this state adiabatically. This 
reduces the effective light-atom interaction dynamics to 
a single decay term 

?f (22|t> <t|p|t) (t|e f - e f c|t> (TIp - Pit) (t| c f c) , (is) 

i.e., the atom is equivalent to a beam splitter which re- 
flects photons out of the cavity at the rate 4g 2 /r if the 
atom is in the state |t) and does not affect the light field 
if the atom is in the state |T). 

Assume the atom to be in the state \j), j G {|,t}- 
Since the input mirror of the cavity may be regarded as 
a beam splitter with high reflectivity, all components of 
the setup transform field operators linearly. For a co- 
herent state input field, the cavity field and the output 
field are hence also coherent states. We may divide the 
time axis into small segments of width r and approximate 
the integrals in (|12[) by sums. The input state is then a 
direct product of single mode coherent states with am- 
plitudes a^tkjyfr and annihilation operators a(tk)y/r, 
where t% = fcr, k = 0, ±1,±2, .... In the following, we 
choose r to be equal to the round trip time of light in 
the cavity, which requires cvi n (t) to vary slowly on that 
time scale. Denoting the coherent state amplitude of the 
cavity field at time t by 7j(£), we use the beam splitter 
transformation 



out 



for the input mirror of the cavity to derive 

Jj(t) = t c y/Ta- m (t) + r c tjJj(t-T), 

V^ a out(*) = r cVra- m {t) - t c tjjj{t - r), 



(16) 



(17) 
(18) 



10 



where r 2 — 1 — t 2 is the reflectivity of the input mirror 
of the cavity, t 2 is the transmissivity of the beam split- 
ter, which models the loss due spontaneous emission, i.e., 
t 2 = 1 — Ag 2 T5j^/T, and a^ t (t) denotes the output field 
from the cavity. We have here included an additional 
phase shift of 7r per round trip in the cavity to ensure the 
input field to be on resonance with the cavity, i.e., the 
second element of the input vector in (|16l) is —tj"f(t — r). 
Taking the limit r — > and t 2 . — > for fixed cavity decay 
rate n := t 2 /r, ((T7J) and HH]) reduce to 

= -(l + 2C<J it )~ 7i (i) + vW*), (19) 



„(i) 



(t) = a in (t) - \/ll^j{t), 



(20) 



where C := 2g 2 /(ftr) is the single atom cooperativity 
parameter. According to the steady state solution of (fl~9|) 



i 



2a in (t) 



1 + 2C5j^ ^/k 



(21) 



we need C 3> 1 to efficiently expel the light field from the 
cavity for j = f. We should also remember the criterion 
for the validity of the adiabatic elimination, which by use 
of (|2Tj) takes the form 

1 » 4g ^in(t)\ 2 ^ K \a in (t)\ 2 ^ 



r 
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i.e., for C>1 the average flux of photons in the input 
beam should not significantly exceed the average number 
of photons emitted spontaneously per unit time from an 
atom, which has an average probability of one half to be 
in the excited state. 
Solving CG 



7j -(i) = V^/ e- (1+2C ^M*-*')/ 2 a in (t')^', (23) 

J — CO 

we obtain the output field o.^ t (t) = a/outM with 

J — oo 

(24) 

Note that f^ t (t) is not necessarily normalized due to the 
possibility of spontaneous emission. The ideal situation 

is /out(*) = /in(*) and f^ t (t) = -/i n (*), and we would 
thus like the norm of 



/oo 
-oo 



(25) 



to be as small as possible. 

For j = j" and C ^> 1, the exponential function in 
(|2"41 is practically zero unless n(t — t') <C 1, and as long 
as /in(i) changes slowly on the time scale (Ck)^ 1 , we 
may take fi n (t) outside the integral to obtain /outW = 
(1 - 2/(1 + 2<7))/m(*) and £ t = 2/(1 + 2C) ~ C" 1 . As 



this result is independent of fi n (t), a natural criterion for 
the optimal choice of input mode function is to minimize 



2 — K 



e- K ^ t,)/2 f* n (t)fin(t')dt'dt 



(26) 

under the constraint f* n (t)fi n (t)dt = 1. We also 
restrict fi n (t) to be zero everywhere outside the time in- 
terval [— T/2,T/2]. Since we would like the double in- 
tegral to be as close to 2 as possible, we should choose 
f* n (t) = fin(t)- A variational calculation then provides 
the optimal solution 



/in(t) 



Acos(uj Q t) for t G [-T/2, T/2] 
otherwise 



where 



.4 



1/2 



T + sm(oj T)/uj 



2lu ( uj T\ 
— tanf — 1=1, 



(27) 



(28) 



(29) 



and € [0, 7r[. 
This solution gives 



E, 



2x 2 



1 



2 cos 2 



2 7 



2^o 

K 



(30) 



For fixed k, wo is a decreasing function of T, and for 
kT > 2tt, w ~ tt/7 1 and 8tt 2 /(kT) 2 . Distortion of 

the input mode function can thus be avoided by choosing 
a sufficiently long input pulse. This may be understood 
by considering the Fourier transform of the input mode 
function. For a very short pulse the frequency distribu- 
tion is very broad and only a small part of the field is ac- 
tually on resonance with the cavity, i.e., most of the field 
is reflected without entering into the cavity. For a very 
long pulse, on the other hand, the frequency distribution 
is very narrow and centered at the resonance frequency of 
the cavity. Ej is plotted in Fig. [S] as a function of T both 
for the optimal input mode function and an input mode 
function, which is constant in the interval [—T/2, T/2]. 
The optimal input mode function is seen to provide a 
significant improvement. In fact, for the constant input 
mode function E± = 4[1 — exp(— kT/2)]/(kT), and hence 
Ei scales only as (kT) -1 for large kT. Finally, we note 
that r = t 2 /n < and kT ^> 1 thus also ensures 

t/T <C 1, which justifies the approximation in Eqs. 1191 
and [201 



B. Single atom cooperativity parameter 

We would like to determine how a single unit of the 
setup transforms the state of the field, when we use the 
full multi-mode description. For this purpose it is simpler 
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FIG. 9: Deviation Ej = 1 + (-l) 6 ^ J f* n (t)f^ t (t)dt of the 
overlap between the output mode function and the input 
mode function from the ideal value when the atom is in the 
state Note that f„{^(t) is defined such that the norm is 
less than unity if a loss occurs during the interaction. The sin- 
gle atom cooperativity parameter is assumed to be C = 10 3 . 
Solid lines are for the optimal input mode function given in 
(|27[) . and dashed lines are for an input mode function, which 
is constant in the interval t £ [— T/2, T/2] and zero otherwise. 
The inset illustrates these functions for kT = 300. For j =f, 
the shape of the output mode function is almost the same as 
the shape of the input mode function as long as kT 3> C _1 , 
but the norm is slightly decreased due to spontaneous emis- 
sion from the atom such that Ef w C" 1 . This result is inde- 
pendent of the actual shape of the input mode function, and 
the solid and dashed lines for Ef are hence indistinguishable 
in the figure. The results for E± are independent of C, be- 
cause the light field does not interact with the atom in this 
case. Nonzero values of E± only occur due to distortion of the 
shape of the mode function. For kT > 1, E± « 8tt 2 /(hT) 2 
for the optimal input mode function and E± m A/(kT) for the 
constant mode function. 



to work in frequency space, and we thus use the definition 
of the Fourier transform f{uS) = J f(t) exp(— iu)t)dt / ^/2/k 
on Eqs. (J23]) and flUJ to obtain 



/ \ _ v 

lj[U> ~ (1 + 2C5jt)«/2 + 

"outM = Kj(w)a- m (w), 



—a- m (oj), 



where 



KM r- 



(1 - 2Cd jt )n/2 - ico 



(31) 
(32) 

(33) 



Assume now that the density operator of the two-beam 
input field to the unit may be written on the form 

Pin = y^y^c» m |{q n (tj)})({a m (a;)}| 

n m 

®\{0n(u)}){{Pm(u)}\, (34) 



where n and m are summed over the same set of num- 
bers and |{a„(w)}) and \{(3 n (uj)}) are continuous coher- 
ent states in frequency space, i.e., 



\{a„(u)}) = exp 



|0) (35) 



and similarly for \{f3 n (uj)}). Note that ([33)1 and (TT21 are 
consistent when a n (ui) and a n (t) as well as a(ui) and a(t) 
are related through a Fourier transform [74j]. After the 
first 50:50 beam splitter, the input state is transformed 
into 



3' = Z! Cnm | + Pn{u))/V2^ 

n 77i 

|{(a m (w)+/3 m H)/\/2} 
{09„(w)-On(w))/^}) 

^{(p m (u)-a m (u)))/V2}\, (36) 



and this is the input state to the two cavities. The initial 
state of the two atoms is 



Pat 



z E E E E m(p\®\j)(i\- (37) 



' ieU,t} ieU.tt pe{l,t} g£{l-t} 



We thus need to know how one of the cavities transforms 
a term like |ai n , n (w))(ai n , m (a;)| (g) \i)(p\. 

We saw in the last subsection that a cavity is equivalent 
to an infinite number of beam splitter operations applied 
to the cavity mode and the input field modes. To take 
the possibility of spontaneous emission into account, we 
also apply a beam splitter operation to the cavity mode 
and a vacuum mode in each time step and subsequently 
trace out the vacuum mode. As the beam splitters are 
unitary operators acting from the left and the right on the 
density operator, the field amplitudes are transformed 
according to (|32|) as before, but when the ket and the bra 
are different, the trace operations lead to a scalar factor. 
Since the reflectivity of the beam splitter modeling the 
loss is 2Ckt8^ for an atom in the state \j), this factor is 

oo 

= Yi ( V 2CKTj myP (kT)5 p ^ I V 2C KTj n ,i(KT)Sif) 
Ck 2 



fc— — oo 



exp 



(1 + 2C) 2 (k/2) 2 + lo- 



■[|a inin (w)| 2 <5 iT 



+ |am,mM| 2 <5 p -t ~ 2a in,n(w)a* nim (w)^^]dw j, (38) 

where J ny i(t) is the amplitude of the cavity field 
corresponding to the input field ot- m;n (t) and the 
atomic state \i) as given in Eq. (|2"3"|) . Altogether, 



12 



l^in.n (w))(ain, m (w)| <& |i)(p| is thus transformed into 

HI® 

|i)(p|. Projecting both atoms onto + |l))/2 and tak- 



ing the final 50:50 beam splitter into account, we obtain 
the output state after one unit 



A>ut=X]X] X! X! X! X! j^ c n m d lp [(a n (uj) + /3 n (u))/V2, (a m (uj) + (3 m (u)) / V2] 

n rn ieU,f} je{-|.,t} p£{l,t} q&U,t} 

x d jq[{l3n(u) - a n (uj))/V2, (P m (u) - a m (uj))/V2} 



-(K p + K q )a m (u) + l(K p ~ K q )(3 m (u) 



^(K p + K q )(3 m (oj) + ^{K p - K q )a m (oj) 



. (39) 



r 



where we have written i"Q = Ki(u>) for brevity. A subse- 
quent phase shifter is easily taken into account by multi- 
plying the coherent state amplitudes by the appropriate 
phase factors. We note that the result has the same form 
as the input state ([34]) if we collect u, i, and j into one 
index and to, p, and q into another index. We can thus 
apply the transformation repeatedly to obtain the output 
state after N + 1 units of the setup. 

If we assume a n (w) and /3 m (uj) to have the same shape 
for all n and to, the output field simplifies to a two-mode 
state for kT > 1 as expected. This may be seen as 
follows. When T is much larger than the width 

of the distribution a n {uj) in frequency space is much 
smaller than k. For the relevant frequencies we thus 
have uj -C fc/2, and in this case the right hand side of 
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FIG. 10: (Color online) Fidelity as defined in ((6} as a function 
of the single atom cooperativity parameter C for different 
values of the expectation value of the number of photons in 
one of the input modes and different numbers of units of the 
setup, assuming kT ^> 1 and pi n = ® l Q )( Q l- The 

dotted lines are the asymptotes for C — > oo. 



p| reduces to -(1 - 2C^ t )/(l + 2C8 jt ), i.e., K 3 is 
independent of frequency, and the amplitudes of all the 
continuous coherent states of the output density operator 
are proportional to a n {ui). 

In Fig. [TTJl w e have used (JSHJ) in the limit kT > 1 
to compute the fidelity © as a function of the single 
atom cooperativity parameter for p- ln = \a)(a\ (g> |cc)(q;|. 
The ideal values for C — > oo are also shown, and the 
deviations are seen to be very small for C larger than 
about 10 3 . Current experiments have demonstrated sin- 
gle atom cooperativity parameters on the order of 10 2 
[l|-y|, and high fidelities can also be achieved for this 
value. When C is finite, there is a relative photon loss 
of 1 — J |/^(i)| 2 di «2/C due to spontaneous emission 
each time the field interacts with a cavity containing an 
atom in the state |f). We note, however, that the setup 
is partially robust against such losses because the next 
unit of the setup removes all components of the state for 
which a single photon has been lost. 



VII. CONCLUSION 

The setup put forward and studied in this article acts 
in many respects like a photon number filter and has 
several attractive applications for quantum technologies. 
Based on the ability to distinguish even and odd pho- 
ton numbers using the interaction of the light field with 
a high finesse optical cavity, photonic two mode input 
states can be projected onto photon-number correlated 
states. Naturally, this protocol is very well suited to de- 
tect losses and can in particular be adapted to purify 
photon-number entangled states in quantum communi- 
cation. We studied deviations from ideal behavior such 
as finite length of the input pulses and limited coupling to 
estimate for which parameters the idealized description is 
valid. The setup can be modified such that it is capable to 
perform a quantum-non-demolition measurement of pho- 
ton numbers in the optical regime. The non-destructive 
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photon counting device completes the versatile toolbox 
provided by the proposed scheme. 
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